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THE INITIAL- VALUE PROBLEM 



The approach to finding exact solutions that has been taken so far has 
involved initially solving the field equations in the interaction region IV 
and then investigating the conditions under which these solutions can 
be considered as the result of collisions of plane waves. In this way, 
resulting solutions are found first and the initial conditions are obtained 
subsequently. In this chapter it is appropriate to return to the original 
problem of specifying the initial data and then attempting to find the 
solution that determines the subsequent development. 

14.1 The initial data 

The problem that is under consideration is the collision of two plane waves 
in a flat background. It is assumed that the two approaching waves are 
both known, and it is required to find the exact solution which describes 
the interaction following the collision. This problem has been formulated 
in earlier chapters. It has been found convenient to divide space-time into 
the four regions as illustrated in Figure 3.1. It has also been demonstrated 
that the metric in the interaction region IV may be taken in the form of 
the Szekeres line element (6.20) which involves the four functions U(u,v), 
V(u, v), W(u,v) and M(u,v) satisfying equations (6.22). One of these 
equations can be integrated to give 

e~ u = f(u)+g(v). (14.1) 

In this approach, it is assumed that the initial conditions are determined 
by the functions f(u), V(u, 0), W(u, 0) and M(it, 0) which represents the 
wave in region II as it reaches the boundary v = 0, and by g(v), V(0, v), 
W(0, v) and M(0, v) which represents the wave in region III as it reaches 
the boundary u = 0. This is now a typical case of the characteristic 
initial-value problem. 

In practice, however, things are a little more complicated since it is 
normal to make use of the transformation (6.7) to put M = in these 
initial regions. In this case the approaching waves are each described by 
three functions either of u or of v satisfying a single equation which is 
either (6.22c) or (6.226). 

The difficulty arises because, as described in previous chapters, it is 
convenient to use / and g, or transformations of them, as coordinates in 
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the interaction region. Accordingly, it is therefore more convenient to use 
the transformation (6.7) in the initial regions to put 

f=\-u n ^{u), g=\-v n *Q{v). (14.2) 

where the constants rt\ > 2 and n 2 > 2 are essential in order to retain the 
continuity properties of the functions / and g across the boundaries v = 
and u = 0. These constants feature prominently in the junction conditions 
as described in Section 7.2. They are determined by the character of the 
wavefronts of the approaching waves. 

In this case, the data specifying the initial waves is now given by 
V(u,0), W(u,0) and M(u,0) on v = 0, and by V(0, v), W(0,v) and 
M(0, v) on u = 0. In addition, since the function M is essentially de- 
termined up to a removable constant for any V and W by equations 
(6.22 b, c), the initial data is effectively described by specifying only the 
metric functions V and W on the boundaries v = and u = 0. It is also 
appropriate to re-express these as functions of / and g. 

The remaining problem is now to integrate the main field equations 
to determine the functions V and W in the interaction region subject to 
their specification on the initial null boundaries. 



14.2 The colinear case 

For the case when the approaching waves have aligned linear polarization 
it is possible to put W = everywhere, and the main field equation for 
V may be expressed in terms of different coordinates in any of the forms 
(9.3), (10.2), (10.11), (10.60) or (10.71). At this point we may take / and 
g as coordinates and consider the equation in the form (10.2) which may 
be rewritten as 

m - v t . + + - o. (".3) 

This equation must now be solved for V(f,g) with initial data defining 
the approaching waves given by V(f, ^), and V(\,g). It is also possible 
to scale the approaching waves such that V = at the point of collision 
u = 0, v = so that V(|, \) = 0. 

Equation (14.3) is an Euler-Poisson-Darboux equation whose solu- 
tion may be expressed as a line integral. As originally pointed out in 
this context by Szekeres (1972) and later repeated by Yurtsever (1988c), 
since (14.3) is a linear hyperbolic equation, it may also be solved explic- 
itly using Riemann's method. According to this method, we consider the 
adjoint equation 
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Figure 14.1 The interaction region IV is represented in f,g coordinates by the 
region inside the triangle shown. The sides g = 1/2 and /= 1/2 are the II— IV 
and III-IV boundaries respectively, and the focusing singularity occurs on the 
line / + g = 0. The solution for V(f,g) may be obtained at any point P by 
integrating round the rectangle PNML. 



where R is a Riemann-Green function satisfying the boundary conditions 
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By using Green's theorem the integral of i?L[V] —VL[R] over the rectangle 
PNML indicated in Figure 14.1 can be related to its line integral around 
the boundary. In this way it can be shown that, for any arbitrary point 
(/o)<7o) within the interaction region, the function V is given by 



V(f ,g )= [ R 

J ML 



V f + 



V 



2(| + /) 



Jmn 



V a + 



V 



dg 
(14.6) 
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provided that R satisfies (14.4) and (14.5) and that V(|, \) = 0. 

A specific Riemann-Green function which satisfies (14.4) and (14.5) 
is given by 



i [ + 9 p A , z(f - fo)(g ~ 9o) i -| i — , 



/o + ^o " 1/Z V (f + g)(fo + g 

where P-1/2 is the Legendre function of order — \. By substituting this 
and the initial expressions for V on the boundaries ML and MN into 
(14.6), an explicit integral for V in the interaction region is obtained. 

This approach of using Riemann's method has been generalized by 
Xanthopoulos (1986) to include the collision of null fluids coupled with 
plane gravitational waves. This particular situation will be discussed 
later in Section 20.2. At this point, we may simply note that the method 
involved is essentially the same. 

Although the method described above does give an explicit integral 
expression for V, in practice it is extremely difficult to evaluate this in- 
tegral for arbitrary initial data. Thus, it does not really provide a viable 
method for generating analytic solutions for any given analytic description 
of the approaching waves. 

Another problem with the method of Riemann is that it is not prac- 
tically possible to generalize it to the case for non-colinear collisions. The 
method only applies to linear hyperbolic equations and, if the polarization 
of the approaching waves is not aligned, the field equations describing the 
resulting interaction are non-linear. In view of this, Hauser and Ernst 
(1989a, 6, 1990) have developed an alternative method which can be gen- 
eralized to the non-colinear situation. 

The method employed by Hauser and Ernst (1989a) is first to find a 
one-parameter family of basic solutions by means of a suitable separation 
of variables and then to express the the final solution as a linear super- 
position of these basic solutions. In fact their approach is based on the 
solution (10.8), which may be re-expressed in the form 

, f{u) A(a) da 



-I 



1/2 \/(J- f(u)y/a + g(v) 

9{v) B(a)da 



+ 



1/2 a/ct + f{u)y/a - gjvj 



for suitable functions A(a) and B(a) which must depend on the initial 
data. 
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An interesting feature of the method is that the functions A(a) and 
B(a) are obtained in terms of the prescribed initial data by solving gen- 
eralized versions of Abel's integral equation for the tautochrone problem 
of classical particle mechanics. Generalized versions of Abel's integral 
equation are required because f'(u) and g'(v ) are zero on the boundaries 
u = and v = respectively. 

Full details of the method have been described by Hauser and Ernst 
(1989a). However, the techniques are still difficult to use in practice for 
arbitrary initial data. 

Hauser and Ernst (19896) have also reformulated the initial-value 
problem for the colinear case as a Hilbert problem in a complex plane in 
two different ways. They have presented solutions of both forms of the 
Hilbert problem and shown that each of these agrees with the solution 
that is obtained by the method referred to above (Hauser and Ernst, 
1989a). 

14.3 The non-colinear case 

Hauser and Ernst (1989a, 6, 1990) have developed the techniques men- 
tioned at the end of the previous section to deal specifically with the case 
when the polarization of the approaching waves is not aligned. 

The initial data have already been described in Section 14.1. It was 
also pointed out that, for the colinear case, it was convenient to rescale 
the approaching waves such that V = at the point of collision. For the 
non-colinear case, Hauser and Ernst have found it convenient to work with 
the complex function Z defined by (11.5) and to use the transformation 
(12.20) which consists of a rotation and rescaling to put Z = 1 at the 
point of collision when u = and v = 0. 

For the colinear case the main field equation (14.3) is linear. However, 
for the non-colinear case this is replaced by two non-linear equations, or by 
the complex Ernst equation. The non-linearity of this equation introduces 
considerable difficulty into any attempt at solving the general problem. 

Hauser and Ernst (1990) have developed a method for considering 
the initial-value problem for general colliding plane gravitational waves. 
They have achieved this by replacing the usual initial-value problem in 
terms of the Ernst equation by an equivalent 2x2 matrix homogeneous 
Hilbert problem in the complex plane. In the case when the polarization of 
the approaching waves is colinear, this approach reduces to the relatively 
simple one-dimensional Hilbert problem that was previously considered 
(Hauser and Ernst, 19896). A detailed description of the method devel- 
oped, however, is beyond the scope of this book and interested readers 
are directed to the original papers. 
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In spite of these advances, there is strong evidence that the general 
solution of the initial-value problem for the non-colinear case is not ex- 
pressible in a finite closed form. In addition, practical methods for solving 
the matrix homogeneous Hilbert problem for particular non-colinear cases 
still need to be developed. 

Clearly this is an area where considerable further work is required. 



